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Abstract 

We present results from a detailed experimental investigation of LaFeAsO, the parent material 
in the series of "FeAs" based oxypnictide superconductors. Upon cooling this material under- 
goes a tetragonal-orthorhombic crystallographic phase transition at ~160 K followed closely by 
an antiferromagnetic ordering near 145 K. Analysis of these phase transitions using temperature 
dependent powder X-ray and neutron diffraction measurements is presented. A magnetic moment 
of ~0.35 hb per iron is derived from Mossbauer spectra in the low temperature phase. Evidence 
of the structural transition is observed at temperatures well above the structural transition (up to 
near 200 K) in the diffraction data as well as the polycrystalline elastic moduli probed by resonant 
ultrasound spectroscopy measurements. The effects of the two phase transitions on the transport 
properties (resistivity, thermal conductivity, Seebeck coefficient, Hall coefficient), heat capacity, 
and magnetization of LaFeAsO are also reported, including a dramatic increase in the magnitude 
of the Hall coefficient below 160 K. The results suggest that the structural distortion leads to a 
localization of carriers on Fe, producing small local magnetic moments which subsequently order 
antiferromagnetically upon further cooling. Evidence of strong electron-phonon interactions in the 
high-temperature tetragonal phase is also observed. 

PACS numbers: 61.05.cp 61.05.fm 61.50.Ks 63.20.kd 65.40.Ba 72.15.Eb 72.15.Jf 74.10. +v 75.25,+z 76.80.+y 
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I. INTRODUCTION 



The family of lanthanide iron oxypnictides crystallizing in the ZrCuSiAs structure type 
has received extensive attention in the recent condensed matter literature. Many of these 
compounds have been known for almost a decade , and superconductivity below ~7 K 
was reported in fluorine doped LaMPO (M = Fe, Ni) in 2006 and 2007 3,3,3- However, the 
discovery of superconductivity at 28 K in fluorine doped LaFeAsO 6J, and the subsequent re- 
ports of transition temperatures greater than 50 K in some of the related rare-earth materials 
has generated great interest and a flurry of recent experimental and theoretical activity, gen- 
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The superconductivity in these materials appears to be unconventional [71, [8|, [29J , and much 
careful work will be required to elucidate the interesting physics in this family of compounds. 
In the arsenic containing materials doping is required to produce the superconducting state, 
and much work is duely focussed on the doping behavior of these compounds. However, 
study of the undoped materials is also important in developing an understanding of the 
superconductivity. 



The undoped material La 



transition near 150 K 



eAsO has been reported to undergo a spin density wave (SDW) 
301 ] . based on specific heat, resistivity, and reflectivity measure- 



ments. A structural phase transition has also been reported in this material at temperatures 



just above the magnetic transition 12 



pressed and superconductivity emerges 



31 



. Upon doping with fluorine the SDW is sup- 



301 ] . Careful characterization of the behavior 



of LaFeAsO and other materials in this family is important in understanding the underlying 
physics responsible for these behaviors. Here we report the results of our experimental in- 
vestigation of LaFeAsO. We present structural analysis through the crystallographic phase 
transition from powder X-ray diffraction (PXRD) and investigation of the magnetic tran- 
sition using neutron powder diffraction (NPD) along with magnetization and Mossbauer 
spectral measurements. Signatures of these phase transitions are observed in ac-calorimetry 
heat capacity and resonant ultrasound spectroscopy (RUS) measurements. The effects of 
the structural and magnetic transitions on the transport properties of LaFeAsO (electrical 
resistivity, magnetoresistance, Hall effect, Seebeck effect, and thermal conductivity) are also 
examined. 
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II. EXPERIMENTAL DETAILS 



A. Sample synthesis 

We tried several methods for synthesizing high quality samples of LaFeAsO. Starting 
materials were from Alfa Aesar and included La (99.9 %), La 2 3 (99.99 %), Fe (99.998 %), 
Fe 2 3 (99.998 %), and As (99.999 %). The binaries LaAs and FeAs were synthesized for use 
as precursors by heating slowly over the course of several days stoichiometric mixtures of 
the elements in sealed, evacuated silica tubes to 950 °C for LaAs and 1050 °C for FeAs, and 
soaking for 12—24 hours. Different synthesis routes will be distinguished by capital letters 
A-D in the discussion below. 

Synthesis route A used a finely ground stoichiometric mixture of LaAs, Fe, and Fe2C>3 
pressed into a pellet, wrapped in Ta foil and heated in a sealed silica tube partially backfilled 
with ultra-high purity Ar and with a small piece of Zr foil, and heated at 1200 °C for 
about one day. Powder X-ray diffraction (PXRD) analysis showed the product to be nearly 
single phase LaFeAsO. The only other phases observed by PXRD were Fe and La2C>3, while 



neutron diffraction experiments detected Fe2As. Rietveld analysis [32] , 1331 . |34| indicates that 
the amount of each of these impurities was less than 5 %. It is likely that the iron containing 
impurities result from the reaction of the Zr foil with As vapor in the tube, reducing FeAs 
(a typical impurity phase when no metal foils are present) to Fe2As and Fe. 

Synthesis route B is identical to A except no metal foils were used. This typically gives 
a product with significant (~10%) La 2 03 and FeAs impurities as determined by neutron 
diffraction and Mossbauer spectroscopy, but no Fe nor Fe 2 As. 

Synthesis route C involves reacting the product of route B with a small excess of La 
powder (a few percent). We found that this helps to further react the La2C>3 and FeAs 
impurities left from route B and can produce PXRD pure material. Occasionally small 
amounts of LaAs and/or La203 are observed in the products from route C. 

Route D used a finely ground stoichiometric mixture of FeAs, La20 3 , and La pressed into 
a pellet, sealed in a silica tube partially backfilled with ultra-high purity Ar, and heated at 
1200°C for 30-36 hours. This method produced purer samples, sometimes with no impurities 
observable by PXRD. 

A single sample prepared by route A was used for neutron diffraction, Hall effect, resistiv- 
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ity, magnetoresistance, thermal conductivity, Seebeck coefficient, and heat capacity measure- 
ments reported below. The presence of elemental Fe precluded magnetic characterization 
of this sample; the material was strongly magnetic at room temperature. Magnetization 
measurements reported below were carried out on a PXRD pure sample obtained via route 
C. Material produced by route B was used for Mossbauer spectroscopy and elastic constant 
measurements. The temperature dependent PXRD data presented here were collected using 
a sample made by route D, and are in good agreement with similar measurements on the 
sample used for transport, NPD, and heat capacity measurements (route A). 



B. Characterization techniques 

Temperature dependent PXRD data were collected using an Anton Parr TTK450 low 
temperature stage on a PANalytical X'pert Pro MPD with an X'celerator position sensitive 
detector and using Cu Ka radiation. The PXRD sample was mixed with copper powder in 
a small amount of vacuum grease to calibrate the sample temperature for each scan using 
the refined lattice constant of Cu, and to increase the thermal contact between the sample 
and the holder. Neutron diffraction measurements were performed at the High Flux Isotope 
Reactor (HFIR) on ~1.5 g samples of LaFeAsO and LaFeAsO0.89F0.11 on the HB1A triple- 
axis spectrometer with horizonal collimations of 48-40-40-68 using a highly filtered beam 
(A/2 ~ 10~ 4 A) with an incident energy of 14.64 meV. Preliminary measurements were also 
performed on the WAND diffractometer and the HB3 triple-axis spectrometer. Data were 
collected over a range of scattering angles from 5° to 130° at several temperatures. 

Transport measurements were performed using a Quantum Design Physical Property 
Measurement System (PPMS). Silver epoxy (Epotek H20E) was used for electrical and 
thermal contacts. Hall effect in fields up to 6 T and electrical resistivity measurements 
were made using platinum wire leads. Gold-coated copper leads were used for thermal 
conductivity and Seebeck coefficient measurements. 

Heat capacity was determined by ac-calorimetry measurements made on a 3.18 mg poly- 



crystalline sample, using 4.5 Hz chopped white light as a heating source [351 ] . The technique 
only yields relative values of the heat capacity, so the results were normalized at 127 K 
to lower resolution specific heat data taken with the PPMS. The data were then corrected 
for the addendum heat capacity (thermocouple wire plus glue), which was ~ 2.5 % of the 
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FIG. 1: (color online) Rietveld refinement results for LaFeAsO at (a) 300 K and (b) 85 K. The 
upper sets of ticks mark the location of Bragg peaks from LaFeAsO, while the lower ticks locate 
Bragg peaks from the Cu internal standard. The inset in (a) shows the tetragonal structure viewed 
along (110) emphasizing the FeAs layer in this material. The grey line outlines the tetragonal 
unit cell. The inset in (b) shows the orthorhombic structure viewed along (0 1), emphasizing 
the distortion which occurs in the ab-plane upon cooling through the tetragonal to orthorhombic 
phase transition. The dashed black lines outline the c-centered orthorhombic unit cell, while the 
dotted grey lines show the primitive unit cell (see text). 



sample's in the temperature region measured. 

The elastic moduli were measured using Resonant Ultrasound Spectroscopy (RUS) and 
a custom made probe in a Quantum Design PPMS cryostat [3^] . Stanford Research synthe- 
sized function generator and Model SR 844 RF Lock-in Amplifier were used to excite the 
sample and collect the data. 

The iron-57 Mossbauer spectra were recorded between 4.2 and 295 K on a constant accel- 
eration spectrometer which utilized a rhodium matrix cobalt-57 source and was calibrated 
at 295 K with a-iron powder. The isomer shifts are relative to room temperature a-iron. 
The accuracy of the temperature is better than ±1 %. 
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III. RESULTS AND DISCUSSION 



A. X-ray diffraction: structural phase transition 

It has been reported that LaFeAsO undergoes a structural phase transition from tetrag- 
onal (P4/nmm, space group number 129) to orthorhombic (Cmme, space group formerly 



known as Cmma, number 67) upon cooling at about 155 K [12l . |31| . Figure [T] shows the 
Rietveld refinement results of powder X-ray diffraction data collected at 300 K and at 85 
K (sample made by synthesis route D). Good fits to the reported structural models are 
observed. The inset of Figure shows the high temperature tetragonal structure viewed 
along the (1 1 0) direction. The layered nature of the structure emphasized in Figure 
persists into the low temperature orthorhombic structure. The inset in Figure [lb shows 
a slab containing one layer of each atom type viewed perpendicular to the layers (along 
the c-axis). Upon cooling through the phase transition the square nets of atoms present 
in the tetragonal structure distort into rectangular nets of Fe and of O, and into centered 
rectangular nets of La and of As. The two Fe-Fe distances in the low temperature phase 
are labeled in Figure [Tb. For comparison, the Fe-Fe distance in the square nets at 300 K is 
2.853 A. 

The structures have been described in detail elsewhere 



12 



31] . However, some confusion 



exists regarding the symmetry of the low temperature structure. It has been reported in 
both primitive monoclinic and c-centered orthorhombic space groups. In the monoclinic 
description the reported lattice is metrically orthorhombic and atoms are very close to 
positions which would give orthorhombic symmetry to the crystal as well. The orthorhombic 
description is probably correct, and no indication of deviation from orthorhombic symmetry 
is observed in our PXRD data. 

To identify the temperature at which the tetragonal-orthorhombic structural transition 
occurs, temperature dependent diffraction data were refined in both models. Figure [2|i 
compares the goodness of fit (x 2 ) obtained from each refinement, and shows that the or- 
thorhombic model gives a better fit at 180 K and below. This suggests that the structural 
transition begins at temperatures significantly higher than previously reported. In addition, 
we compare in Figure [2^ the width of two Bragg peaks, one which is split by the structural 
distortion (110) and one which is not (2 0). We note that indices refer to the tetragonal 
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structure. The results show that while the (2 0) peak width is unchanged upon cooling, 
the (1 1 0) peak begins to broaden significantly at about 180 K. 

Figure [2b shows the temperature dependence of the refined lattice constants. Deter- 
mining the temperature at which to switch between the Cmme and P4/nmm models is not 
straightforward. We chose to show data in the Cmme model when the difference between the 
refined values of the a and b axes is greater than the error associated with these parameters. 
Based on this criterion, the Cmme model was used at 200 K and below. Above this tem- 
perature refinement in this model is in fact unstable and P4/nmm was used. Inspection of 
the refined values of the lattice constants a and b reveals upon cooling a gradual divergence 
below 200 K followed by a more rapid divergence at lower temperatures. The temperature 
dependence of the lattice constant c is shown in the inset of Figure [2b, and undergoes a 
change in slope through the structural transition. 

The above observations suggest that the structural distortion in LaFeAsO occurs over a 
wide temperature range, but includes a "sharp" anomaly as well. To derive a transition tem- 
perature Tt-o from the structural data, we plot the difference between the lattice constants 
a and b in Figure [2b. These data show an abrupt slope change at about 160 K which we 
identify as Tt-o- We note however that the onset temperature of this distortion (180—200 
K) is significantly higher than Tt-o- We believe this extension of the lattice distortion to 
temperatures well above the "transition temperature" is observable in the data shown in 



previous reports which interpreted the results as a single sharp transition jl2|, [31j. Per- 
haps the two-dimensional nature of the crystal and electronic structures leads to enhanced 
fluctuations above the transition temperature. 

B. Neutron diffraction: magnetic phase transition 

Neutron diffraction results reveal a magnetic transition somewhat below the temperature 
at which the structural transition occurs. Figure [3b shows neutron diffraction data for both 



LaFeAsO (synthesis route A) and superconducting LaFeAsOo.sgFo.n |s| at 70 and 200 K 
from 32 to 36° in scattering angle. These data indicate the presence of additional scattering 
at 33.7° at 70 K that is not present at 200 K. Figure [3b shows the temperature dependence 
of the scattering at 33.7°. These data show that the additional low temperature intensity 
in LaFeAsO appears below 145 K. Polarized measurements (not shown) indicate that this 
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FIG. 2: (color online) Results of temperature dependent PXRD analysis of LaFeAsO showing the 

continuous nature of the crystallographic phase transition, (a) Solid symbols show x 2 from Rietveld 

refinement using the orthorhombic low temperature model and tetragonal high temperature model. 

Open symbols show the normalized width of two Bragg peaks labeled by their tetragonal indices. 

The (110) peak is expected to be split into two peaks by the orthorhombic distortion, while 

the (2 0) is not. (b) Unit cell parameters determined by Rietveld refinement. For comparison 

the a lattice parameter in the tetragonal structure has been multiplied by y/2. Error bars are 

not shown on some data points for which the error is smaller than the size of the data point 

symbol, (c) The difference between b and a in the orthorhombic model, linearly coupled to the 

order parameter for the structural phase transition. We define the transition temperature for the 

tetragonal to orthorhombic transition as the kink near 160 K, but note that the structure is evolving 
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continuously both above and below this temperature. 
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FIG. 3: (color online) (a) Neutron diffraction data for LaFeAsO and LaFeAsO0.89F0.11 at 70 and 
200 K. Arrows indicate the angular position of the (1 0) and (0.5 0.5 1.5) wavevectors. For clarity, 
the data for LaFeAsO has been displaced by 150 counts, (b) The temperature dependence of the 
(0.5 0.5 1.5) and (0.5 0.5 0.5) (inset) positions. The wavevectors are labeled using the tetragonal 
setting. The lines are a guide to the eye. 

scattering is magnetic in nature. The same angular range was explored in LaFeAsO0.89F0.11- 
Although the counting time was doubled to 4 minutes per point (normalized in Fig. [3^ to 2 
minutes per point for the purposes of comparison) there is no discernable difference between 
the data at 70 and 200 K. Full refinements of the neutron data indicate similar impurity 
phases in both samples and, hence, the extra intensity is intrinsic to LaFeAsO. For reference, 
the strong rise in intensity below 32.5° corresponds to the (0 2) structural Bragg peak. 

The arrows in Fig. [3]^a) show the position of two plausible wavevectors (labeled in the 
tetragonal setting) to describe the additional intensity at 33.7° in LaFeAsO. The position 
of (1 0) is significantly far from the observed peak in the intensity to rule it out as the 
wavevector. However, the data is strongly consistent with a wavevector of (0.5 0.5 1.5). 
To check the indexing of this peak, we searched for extra intensity at (0.5 0.5 0.5). The 
intensity observed at this location is very weak but careful measurement of the tempera- 
ture dependence shown in the inset of Fig. [3]^b) indicates an increase below 145 K in a 
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FIG. 4: (color online) Results of resonant ultrasound spectroscopy studies showing the effects of 
the crystallographic phase transition on the polycrystalline elastic moduli en and C44 of LaFeAsO. 
A gradual softening begins above 200 K, well above the structural transition, and extends down to 
140 K. A relatively sharp peak in the derivative (inset) occurs near Tt~o- 

manner consistent with the 33.7° reflection. The observation of these two reflections pro- 
vides strong evidence that the true ordering wavevector in LaFeAsO is (0.5 0.5 0.5). This 
indicates a doubling of the conventional unit cell along both a-axes and along the c-axis, 
and is also consistent with the y/2 x y/2 unit cell suggested by Dong et al. 30[], and is in 
good agreement with the magnetic structure derived from other recent neutron diffraction 
experiments on LaFeAsO 12j. Evidence of this antiferromagnetic ordering is also observed 
in the magnetization data discussed below. 



C. Elastic constants 



Resonant ultrasound spectroscopy (RUS) measurements were performed to investigate 
the elastic properties of LaFeAsO. This technique can be a very sensitive and powerful 
probe for studying phase transitions, in particular structural transitions in which strain is 
the order parameter. Figure H] shows the temperature dependence of the polycrystalline 
elastic moduli en and C44. A remarkably gradual softening is seen to extend up to above 200 
K. Data are shown for a sample made by synthesis route B; similar behavior was observed in 
material from synthesis route A. The temperature derivative of the elastic moduli are shown 
in the inset of Figure HI The sharp cusp in dc i: ,/dT occurs at 156 K and is identified as 
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FIG. 5: (color online) The measured heat capacity of LaFeAsO (in units of R per formula unit) 
which displays broad anomalies at both phase transitions. Inset (a) shows the baseline subtracted 
data which clearly shows two broad peaks centered near the magnetic and structural transition 
temperatures. Inset (b) shows the integrated entropy associated with the two anomalies. 



the structural transition temperature determined by this experimental probe. This agrees 
well with Tt-o determined by the structural analysis presented above. The extension of 
the elastic softening to well above Tt-o is consistent with the PXRD data presented above 
and is further evidence of the gradual nature of the crystallographic phase change in this 
material. We also note that no significant thermal hysteresis or magnetic field effect on the 
elastic properties was observed. 



D. Heat capacity 

The measured heat capacity of LaFeAsO (synthesis route A) in the vicinity of the struc- 
tural and magnetic transitions is shown in Figure [5j Two broad overlapping anomalies are 
observed. We note that two heat capacity anomalies were also observed in sample used 
for PXRD analysis, produced by route D. Baseline subtracted data are shown in Figure [5] 
inset (a). The background was estimated by a polynomial fit to the data above and below 
the transition region. The subtracted data clearly show two peaks, one associated with the 
structural transition centered at 155 K, and one due to the magnetic transition centered at 
143 K. The coincidence of these two peaks with the structural and magnetic phase tran- 
sition temperatures strongly suggest that these are indeed separate anomalies, and not a 
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FIG. 6: (color online) Effects of the phase transitions on the electrical transport of LaFeAsO. 
(a) The temperature dependence of the electrical resistivity with no applied magnetic field and 
with an applied field of 8 T. The inset in (a) shows the magnetoresistance calculated from the 
resistivity data, (b) The temperature derivative of the measured resistivity on cooling and warming 
illustrating the absence of thermal hysteresis. The effect of the structural transition at 160 K is 
shown. The peak in dp/dT is near the magnetic transition temperature at 143 K. 

single transition smeared by, for example, inhomogeneities. No hysteresis was observed in 
repeated measurements, suggesting that the phase transitions are second order or perhaps 
weakly first order. The entropy change determined by integration of the subtracted heat 
capacity data is shown in Figure [5] inset (b). A total entropy of 0.032 R (0.27 J K _1 mol 
FU _1 ) is determined by the integration. If the total change in entropy AS is considered 
purely electronic (AS = TqA'j) then A7 ~1.7 mJ mol -1 K -2 . 



E. Transport properties 

The transport properties reported here were all measured on a single sample produced 
using synthesis route A. Figure [6^ shows the measured electrical resistivity of LaFeAsO, 
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which agrees well with previous reports [37|. At room temperature p has a value of 4 
mQcm and decreases upon cooling. This is typical of a low carrier concentration metal 
or heavily doped semiconductor. The upturn in p on cooling below about 200 K is likely 
due to increased charge carrier scattering by lattice fluctuations related to the onset of the 
impending structural transition. The electrical resistivity reaches a local maximum at 165 
K and drops rapidly below this temperature. 

The temperature derivative of the measured resistivity data is shown in Figure [6b. Below 
about 160 K a sharp increase in dp /dT is observed. The maximum occurs at 143 K, 
near the magnetic transition temperature. No thermal hysteresis is observed in these data, 
suggesting that neither of the phase transitions are strongly first order. Since the resistivity 
anomaly begins near Tt-o — 160 K and this is the temperature below which significant 
magnetoresistance appears (Figure [6^ inset), we believe that structural phase transition in 
LaFeAsO is primarily responsible for the dramatic changes in transport properties. The 
magnetic transition occurs near the peak in dp/dT, and may be responsible for the weak 
anomaly observed at 142 K in the magnetoresistance (Figure Eh inset). 

Results of Hall effect measurements are shown in Figure [7J Near room temperature the 
Hall coefficient is negative and nearly temperature independent, indicating conduction by 
electrons with an inferred concentration n = 3 x 10 21 cm -3 . We note that the presence 
of multiple bands at the Fermi level complicates the interpretation of the Hall coefficient 
data. In a simple model with one electron band and one hole band the inferred value of n 
gives an upper bound on the concentration of the dominant charge carriers. Upon cooling 
below 160 K the Hall coefficient decreases rapidly, and n drops by an order of magnitude 
between 160 and 100 K (Figure [TJi, inset). This suggests that many of the charge carriers 
present in the high temperature phase are localized at Tt~o due to the structural transition. 
This is possibly related to the local moment formation and subsequent magnetic ordering 
observed below Tt-o- It is interesting that this order of magnitude decrease in n coincides 
with a factor of two decrease in p. This suggests a large change in carrier mobility. The 
Hall mobility fin = RhP" 1 is plotted in Figure [7b. The mobility increases by a factor 
of 300 between room temperature and 5 K, and reaches the remarkably high value of 150 
cm 2 V~ l s^ 1 in this polycrystalline material. 

The measured thermal conductivity k and Seebeck coefficient S of LaFeAsO are shown 
in Fig. [8] The thermal conductivity increases abruptly below about 155 K, but otherwise 
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FIG. 7: (color online) Results of Hall effect measurements on LaFeAsO showing the remarkable 
decrease in inferred carrier concentration and high mobility at low temperatures, (a) The measured 
Hall coefficient and inferred carrier concentration (inset), (b) The Hall mobility calculated from 
the carrier concentration in (a) and the electrical resistivity (Figure E]) . 

follows the behavior of typical crystalline materials. The increase in k cannot be attributed 
to the change in electronic thermal conduction as the thermal conductivity in this low 
carrier concentration material is dominated by phonons. Thus the observed behavior must 
be attributed to an increase in the thermal conductivity of the lattice, and related to the 
tetragonal-orthorhombic crystallographic transition that occurs at 160 K (Figure[H] inset a). 
This could be due to a decrease in electron-phonon scattering below the transition, which 
would suggest strong coupling between the charge carriers and the lattice vibrations in the 
tetragonal phase of LaFeAsO through bond length fluctuations. This is also consistent 
with the rapid increase in the carrier mobility presented above (Figure [7}d). An increase 
in phonon thermal conductivity could also arise from the freezing in of phonon-scattering 
lattice fluctuations upon cooling through T T _ . 

The Seebeck coefficient (Fig. [8] inset b) is negative over the entire temperature range, 
indicating that electrons dominate the electrical conduction. This is consistent with the 
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FIG. 8: (color online) Thermal and thermoelectric transport properties of LaFeAsO. The thermal 
conductivity measured in zero applied magnetic field and in a field of 8 T, showing an abrupt 
change in slope (inset a) as the structure transforms from tetragonal to orthorhombic upon cooling. 
Inset (b) shows the dramatic changes in the Seebeck coefficient that occur below the transition 
temperatures. 

negative Hall coefficient (Figure [7]). The Seebeck coefficient is moderately high in this 
material, but lower by about a factor of two than in superconducting LaFeAsO0.89F0.11 
above Tq [8]. A remarkable decrease in 15*1 is observed below about 155 K. It is unusual to 
see a sharp drop in \S\ coinciding with a sharp drop in carrier concentration. Boltzmann 



transport theory predicts for the free electron model 



38 



S(T) 



3 lei 



k B T 



+ 



dr 



11 



(1) 



The observed decrease in both \S\ and n through the transition indicates that the second 
term in Eq. [T] is dominant in this temperature regime. This suggests that the charge 
carrier scattering mechanism is changed significantly as the material passes through the 
phase transition region, and presents further evidence for the reduction of electron-phonon 
interactions in the orthorhombic phase which was suggested by the above analysis of carrier 
mobility and thermal conductivity. This is evidence of strong electron-phonon coupling in 
LaFeAsO. 
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FIG. 9: (color online) The temperature dependence of M/H (per mole of formula units) for a PXRD 
pure polycrystalline LaFeAsO sample. The inset shows M(H) data at four different temperatures. 

F. Magnetic properties 

Results of magnetization measurements on a PXRD pure polycrystalline sample of 
LaFeAsO (synthesis route C) are shown in Figure [91 Measured M(H) curves (Fig. [9] in- 
set) show a paramagnetic response up to H = 6 T at 100 K, 150 K, and 200 K. The upturn 
in M/H at low temperatures and the nonlinear behavior of M(H) at 1.8 K are likely due 
to small amounts of impurities. Other than the Curie tail at low temperatures, the mag- 
netization shows little temperature dependence above and below the transition region. A 
decrease in M/H is observed through the phase transitions, beginning near 153 K. This has 
been identified as a signature of spin density wave formation. Localization of charge carri- 
ers which occurs at T^-o would also result in a decrease in M/H through the reduction of 
the Pauli paramagnetism; however, one would expect Curie- Weiss behavior below this tem- 



perature i 
elsewhere 



that were the complete story. Neutron diffraction results presented above and 



12| suggest an antiferromagnetic ordering develops near the temperature at which 



M/H decreases in this material. Carrier localization followed closely by antiferromagnetic 
ordering could be responsible for the behavior observed near 153 K in Fig. [91 The magnetic 
behavior of this anisotropic material is certainly complex and not yet well characterized or 
well understood. Future studies of single crystals with well controlled stoichiometry will be 
of great importance to the understanding of the magnetic nature of LaFeAsO. 

The Mossbauer spectra of LaFeAsO produced by synthesis route B are shown in Fig. [TU1 
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The spectra were modeled using a mixture of LaFeAsO and FeAs (vide infra). Between 4.2 
and 75 K, the LaFeAsO Mossbauer spectrum is a simple magnetic sextet. At 4.2 K, the 
isomer shift and quadrupole shift of LaFeAsO are 0.576(5) and -0.031(1) mm/s, respectively, 
indicating the low-spin nature of the iron(II). The hyperfine field is 5.19(1) T and the usual 



conversion factor of 15 T per 1 /ig yields an estimated Fe moment of ca. 0.35 
hyperfine parameters are in very good agreement with those reported in Ref. 



Mb 



39 



40 



ncsc 



|41|. 
42]. 



A recent theoretical study has also derived 0.25-0.35 in the low temperature phase 
Note that in Ref. [4^] the LaFeAs sample was oxidized and that the spectra are actually 
those of LaFeAsO. The small, but non-zero, quadrupole shift indicates that there is a small 
lattice contribution to the iron(II) quadrupole interaction, as expected for low-spin iron(II) 
in a slightly distorted tetrahedral environment. The observed quadrupole shift is constant 
up to 140 K. At 150 and 295 K the quadrupole splitting is zero within the error bar. The 
introduction of a small, temperature independent, texture for the LaFeAsO phase, with 
intensity ratios of 3:2.3:1:1:2.3:3 yielded a significant reduction in % 2 and this texture was 
used for the fits in Fig. [101 The temperature dependence of the isomer shift is in very good 
agreement with the Debye model for the second order Doppler shift and yielded a Mossbauer 
temperature of 294(20) K, see lower inset in Fig. [TUJ The magnetic moment is essentially 
constant up to 75 K, and the gradual reduction of the hyperfine field above 125 K indicates 
a transition to a paramagnetic phase close 150 K. Between 125 and 150 K, the LaFeAsO 
spectrum is more complex, as in Refs. 37|, [39J], and was modeled herein as a superposition 
of a sextet and a singlet, and the weighted average hyperfine field is shown in the upper 
inset of Fig. [TOJ The observed spectra between 125 and 150 K could be explained either 
by a small compositional variation and a smearing of or by a hyperfine field distribution 



37 



39] or from spin 



resulting from a incommensurate or commensurate spin density wave 
glass behavior in this temperature range. At 75 K and below, however, the simple sextet 
spectrum is indicative of no such spin density wave or spin glass behavior. Finally, we 
observe a marginal increase in linewidth on cooling from 295 to 150 K, this increase could 
be related to the tetragonal to orthorhombic structural distortion. 

Because the Mossbauer spectra of LaFeAsO and FeAs 43|, |44j are located in the same 
velocity region, detecting FeAs impurities requires careful comparison of the spectra near the 
77 K Neel temperature of FeAs. It is clear that the transition of FeAs from the helimagnetic 
to paramagnetic phase is responsible for the modification of the spectra between 65 and 75 
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K. The FeAs contribution above and below 70 K was found to be consistent with and then 
constrained to the hyperfine parameters for the corresponding temperatures of Ref. 
Ref. [43j, respectively, and resulting fits are of very good quality. The FeAs contribution 



44j and 



corresponds to ca. 10 % weight of the sample, a quantity that is essentially temperature 
independent, indicating a similar temperature dependence of the recoil free fraction in FeAs 
and LaFeAsO and indeed the ca. 285 K Mossbauer temperature obtained from the second 



order Doppler shift in FeAs, Ref. 



44J and the line in the lower inset in Fig. [TOj is very 



close to 294 K observed herein for LaFeAsO. The only free fit parameters for the FeAs phase 
are the linewidth, and the isomer shift and average hyperfine field at 25, 50, and 65 K. We 



3,Q 



arises not 



believe that the larger intensity of the 2:5 lines in the spectra of Refs. [37J, 
from a possible texturing of the sample but from an underlying FeAs impurity subspectrum. 
The presence of this FeAs imp urity , if not properly modeled will lead to a smaller effective 
hyperfine field, as seen in Ref. Bother, it is likely that the third component observed 
in /iSR measurements below 70 K [37J is also related to the FeAs impurity that is clearly 
visible in the 78 K Mossbauer spectrum on the same sample. 



IV. CONCLUSIONS 



At high temperatures LaFeAsO is a low carrier concentration metal with conduction dom- 
inated by electrons and with no local magnetic moment. This is supported by measurements 
of electrical resistivity (Figure [6]), Hall coefficient and carrier concentration (Figure [7]), See- 
beck coefficient (Figured]), and magnetization (Figure [9]). Upon cooling a crystallographic 
phase transition occurs. We propose that the structural transformation occurs continuously 
below about 200 K based on structure refinements (Figure [2]) and elastic response (Figure 
Hj), and that the kink in the order parameter b-a at 160 K (Figure [2b ) indicates the point 
at which sufficient distortion has occurred to result in carrier localization and local moment 
formation on the Fe atoms. This is supported by the small entropy associated with the 
transition at Tt-o (Figure [5]). It is unclear what drives this structural transition. One 
candidate is a band Jahn- Teller effect j^j], in which energy is gained by splitting the sharp 
peak in the density of states near the Fermi level [h]] by distorting the crystal structure. 
This could result in localization of some conduction electrons and the development of a local 
magnetic moment on Fe. Careful theoretical analysis of the effects of structural distortion 
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FIG. 10: (color online) The Mossbauer spectra of LaFeAsO. The total data and fits, as well as 
the data with the FeAs impurity contribution substracted and the LaFeAsO phase fit are shown. 
Upper inset: the average hyperfine field determined in the LaFeAsO phase. Lower inset: The 
temperature dependence of the isomer shift of LaFeAsO, with a Debye model fit, and, straight 
line, the isomer shift in FeAs from Ref. [l4]. 
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on the band structure of LaFeAsO could help resolve this issue. Upon further cooling the 
local moments order near 145 K (Figure [3]) . Strong electron-phonon coupling exists in the 
high temperature tetragonal phase, as evidenced by the behavior of the mobility (Figure[7b), 
thermal conductivity, and Seebeck coefficient (Figured]) through the phase transition region. 
Upon doping with, for example, fluorine, the structural phase transformation and associated 
magnetic ordering is suppressed (Figure [3]), allowing the strong electron-phonon coupling 
present in the tetragonal phase to extend down to low temperature, and superconductivity 
emerges. It is believed that the phase transition is suppressed electronically by electron 
or hole doping with substitutional F or Sr, or with oxygen vacancies. However, disorder 
may also suppress the structural transition. Perhaps new superconductors can be discov- 
ered by investigating isoelectronic substitutions which suppress Tt-o i n the ZrCuSiAs-type 
rare-earth iron oxyarsenides and related ThCr2Si2-type alkaline-earth iron arsenides. 

We acknowledge useful discussions with D. J. Singh, O. Garlea, S. Nagler, A. Castro- 
Neto, and I. I. Mazin, and technical assistance by J. Zarestky and J. Fernandez-Baca. 
Research sponsored by the Division of Materials Sciences and Engineering, Office of Basic 
Energy Sciences. Part of this research performed by Eugene P. Wigner Fellows at ORNL, 
managed by UT-Battelle, LLC, for the U.S. DOE under Contract DE-AC05-00OR22725. 
A portion of this research at ORNL's Center for Nanophase Materials Sciences and High 
Flux Isotope Reactor was sponsored by the Scientific User FacilitiesDivision, Office of Basic 
Energy Sciences, DOE. A portion of this work was supported by NSF grant DMR-0400938. 
F.G. acknowledges with thanks the financial support of the Fonds National de la Recherche 
Scientifique, Belgium, through grants 9.456595 and 1.5.064.05. 



[1] P. Quebe, L. J. Terbiie, and W. Jeitschko, J. Alloys Compd. 302, 70 (2000). 
[2] B. Zimmer, W. Jeitschko, J. H. Albering, R. Glaum, and M. Reehuis, J. Alloys Compd. 229, 
238 (1995). 

[3] Y. Kamihara, H. Hiramatsu, M. Hirano, R. Kawamura, H. Yanagi, T. Kamiya, and H. Hosono, 

J. Am. Chem. Soc. 128, 10012 (2006). 
[4] C. Y. Liang, R. C. Che, H. X. Yang, H. F. Tian, R. J. Xiao, J. B. Lu, R. Li, and J. Q. Li, 

Supercond. Sci. Lechnol. 20, 687 (2007). 



21 



[5] T. Watanabe, H. Yanagi, T. Kamiya, Y. Kamihara, H. Hiramatsu, M. Hirano, and H. Hosono, 

Inorg. Chem. 46, 7719 (2007). 
[6] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. Am. Chem. Soc. 130, 3296 (2008). 
[7] F. Hunte, J. Jaroszynski, A. Gurevich, D. C. Larbalestier, R. Jin, A. S. Sefat, M. A. McGuire, 

B. C. Sales, D. K. Christen, and D. Mandrus, Nature 453, 903 (2008). 
[8] A. S. Sefat, M. A. McGuire, B. C. Sales, R. Jin, J. Y. Howe, and D. Mandrus, Phys. Rev. B 

77, 174503 (2008). 

[9] X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen, and D. F. Fang, Nature 453, 761 (2008). 
[10] D. J. Singh and M. H. Du, Phys. Rev. Lett. 100, 237003 (2008). 
[11] Y. Sun, M. Guidry, and C. Wu, Chinese Science Bulleitn 53, 1617 (2008). 
[12] C. de la Cruz, Q. Huang, J. W. Lynn, J. Li, W. Ratcliff II, J. L. Zarestky, H. A. Mook, G. F. 

Chen, J. L. Luo, N. L. Wang, et al., Nature 453, 899 (2008). 
[13] D. Wu, Z. Li, S. Su, J. Dong, G. Li, W.-Z. Hu, P. Zheng, G.-F. Chen, J.-L. Luo, and N.-L. 

Wang, Science in China Series G 51, 715 (2008). 
[14] P. Cheng, L. Fang, H. Yang, X.-Y. Zhu, G. Mu, H.-Q. Luo, Z.-S. Wang, and H.-H. Wen, 

Science in China Series G 51, 719 (2008). 
[15] S. Ishibashi, K. Terakura, and H. Hosono, Journal of the Physical Society of Japan 77, 053709 

(2008). 

[16] C. Day, Phys. Today 61, 11 (2008). 

[17] H. Takahashi, K. Igawa, K. Arii, Y. Kamihara, M. Hirano, and H. Hosono, Nature 453, 376 
(2008). 

[18] Q. Han, Y. Chen, and Z. D. Wang, EPL 82, 37007 (2008). 
[19] K. Haule, J. H. Shim, , and G. Kotliar, Phys. Rev. Lett. 100, 226402 (2008). 
[20] H.-H. Wen, G. Mu, L. Fang, H. Yang, and X. Zhu, EPL 82, 17009 (2008). 
[21] E. Hand, Nature 452, 922 (2008). 

[22] H. Kito, H. Eisaki, and A. Iyo, J. Phys. Soc. Japan 77, 063707 (2008). 

[23] T. Sato, S. Souma, K. Nakayama, K. Terashima, K. Sugawara, T. Takahashi, Y. Kamihara, 

M. Hirano, and H. Hosono, J. Phys. Soc. Japan 77, 063708 (2008). 
[24] S. Ishibashi, K. Terakura, and H. Hosono, J. Phys. Soc. Japan 77, 053709 (2008). 
[25] T. Y. Chen, Z. Tesanovic, R. H. Liu, X. H. Chen, and C. L. Chien, Nature, 

doi: 10. 1038 /nature07081 . 

22 



[26] Z.-A. Ren, J. Yang, W. Lu, W. Yi, X.-L. Shen, Z.-C. Li, G.-C. Che, X.-L. Dong, L.-L. Sun, 

F. Zhou, et al., EPL 82, 57002 (2008). 
[27] Z.-A. Ren, J. Yang, W. Lu, W. Yi, G.-C. Che, X.-L. Dong, L.-L. Sun, and Z.-X. Zhao, EPL 

83, 17002 (2008). 

[28] G. Xu, W. Ming, Y. Yao, X. Dai, S.-C. Zhang, and Z. Fang, EPL 82, 67002 (2008). 

[29] I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, Phys. Rev. Lett. 101, 057003 (2008). 

[30] J. Dong, H. J. Zhang, G. Xu, Z. Li, G. Li, W. Z. Hu, D. Wu, G. F. Chen, X. Dai, J. L. Luo, 

et al., cond-mat/0803.3426. 
[31] T. Nomura, S. W. Kim, Y. Kamihara, M. Hirano, P. V. Sushko, K. Kato, M. Takata, A. L. 

Shluger, and H. Hosono, cond-mat/0804.3569. 
[32] J. Rodriguez-Carvajal, FullProf Suite 2005, version 3.30, June 2005, ILL. 
[33] B. H. Toby, J. Appl. Cryst. 34, 210 (2001). 

[34] A. C. Larson and R. B. V. Dreele, General Structural Analysis System (GSAS), Los Alamos 

National Laboratory Report LAUR 86-748 (2004). 
[35] M. Chung, E. Figueroa, Y.-K. Kuo, Y. Wang, J. W. Brill, T. Burgin, and L. K. Montgomery, 

Phys. Rev. B 48, 9256 (1993). 
[36] A. Migliori and J. Maynard, Rev. Sci. Instrum. 76, 121301 (2005). 

[37] H.-H. Klauss, H. Luetkens, R. Klingeler, C. Hess, F. Litterst, M. Kraken, M. M. Korshunov, 

I. Eremin, S.-L. Drechsler, R. Khasanov, et al., Phys. Rev. Lett. 101, 077005 (2008). 
[38] J. McCarten, S. E. Brown, C. L. Seaman, and M. B. Maple, Phys. Rev. B 49, 6400 (1994). 
[39] S. Kitao, Y. Kobayashi, S. Higashitaniguchi, M. Saito, Y. Kamihara, M. Hirano, T. Mitsui, 

H. Hosono, and M. Seto, cond-mat/0805.0041. 
[40] H. Ramus, E. Morsen, B. Mosel, W. Miiller-Warmuth, W. Jeitschko, L. Terbiichte, and 

T. Vomhof, J. Phys. Chem. Solids 54, 135 (1993). 
[41] I. Nowik, I. Felner, V. Awana, A. Vajpayee, and H. Kishan, J. Phys. Condens. Matter 20, 

292201 (2008). 

[42] J. Wu, P. Phillips, and A. H. Castro Neto, cond-mat/0805.2167. 

[43] L. Haggstrom, A. Gustavsson-Seidel, and H. Fjellvag, Europhys. Lett. 9, 87 (1989). 

[44] S. Kulshreshtha and P. Raj, J. Phys. F: Metal Phys. 9, 2253 (1979). 

[45] R. N. Bhatt, Phys. Rev. B 16, 1915 (1977). 



23 



